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Abstract 

It is shown that the percolation model of hydrogen-bonded crys- 
tals, which is a 6-vertex model with bond defects, is completely equiv- 
alent with an 8-vertex model in an external electric field. Using this 
equivalence we solve exactly a particular 6-vertex model with bond 
defects. The general solution for the Bethe-like lattice is also ana- 
lyzed. 



1 Introduction 

The 6-vertex model on a square lattice describes hydrogen-bonded crystals 
in two dimensions. Historically, it was Slater who first considered the 
evaluation of the residual entropy of ice, a hydrogen-bonded crystal, under 
the assumptions that (i) there is one hydrogen atom on each lattice edge, 
and (ii) there are always two hydrogen atoms near, and away from, each 
lattice site (the ice rule). Under these assumptions there are 6 possible 



hydrogen configurations at each site, and one is led to a 6-vertex model. The 
exact residual entropy of the "square" ice, i.e., ice on the square lattice, was 
obtained by Lieb which gives rise to a numerical number surprisingly 
close to the experimental residual entropy of real ice (in three dimensions). 
The 6-vertex model is therefore an accurate description of hydrogen-bonded 
crystals. 

In real hydrogen-bonded crystals, however, there exist bonding defects 
[|]. One way through which bond defects can occur is caused by the double- 
well potential seen by hydrogen atoms between two lattice sites. When two 
hydrogens occupy the two off-center potential wells along a given lattice 
edge, the assumption (i) above is broken, albeit the ice rule (ii) is still intact. 
This leads to the fact that one, two or zero hydrogen atoms can be present 
on a lattice edge. Indeed, one of us has considered this possibility in 
a percolation model of supercooled water. The same model has later been 
considered by Attard and Batchelor who analyzed it using series analyses. 
More recently, Attard reformulated the problem as a 14- vertex model with 
Bjerrum bond defects, and analyzed the 14-vertex model using an indepen- 
dent bond approximation. Here, using a somewhat different mapping, we 
establish the exact equivalence of the 6-vertex model with bond defects with 
an 8- vertex model in an external field. As a result, we are able to analyse 
the exact solution in a particular parameter subspace. We also discuss the 
general solution of the 6-vertex model with bond defects on the Bethe lattice. 

2 Equivalence with an 8-vertex model 

Consider a square lattice £ of sites under periodic boundary conditions so 
that there are 2N lattice edges. The lattice is hydrogen-bonded with defects 
such that there can be one, two or zero hydrogen atoms on each lattice edge. 
As the hydrogen atoms are placed off-center on the edges, we place two Ising 
spins 0", a' on each lattice edge such that a = 1 denotes that the site is 
occupied by a hydrogen and a = —1 the site is empty. However, the ice rule 
dictates that the sum of the 4 Ising spins ai, (J2, as, surrounding a square 
lattice sites must vanish. There are altogether six ice-rule configurations as 
shown in Fig. 1. 

More generally, we consider the ice-rule model with weights 

uJi = UJ2 = a, uo^ = 00^ = 6, 00^ = ujq = c, (1) 
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Figure 1: The six ice-rule configurations. Each sohd circle denotes an Ising 
spin a = 1 and each open circle a a = — 1. 

where Ui is the weight of the ith configuration shown in Fig. 1. Denote the 
weights (|lD by ci;(cri, cr2, o"3, (T4) where the subscripts are indexed as shown 
and ui = uj{l, —1, —1, 1), etc. The weight uj{ai, (T2, cts, a^) satisfies the "spin 
reversal" symmetry 

uj{ai, (72, (T3, (T4) = -0-2, -0-3, -(^4) (2) 

and vanishes except for the six weights given in (|l]). To each lattice edge 
containing two spins a and a', introduce an edge factor E{a, a') to reflect 
the effect of bond defects. Then, the partition function of interest is 

Z = n ^^(cri,0-2,0-3,Cr4) n ^(^'^')- (3) 

o-=±l vertices edges 

In the 6-vertex model without bond defects we have E{a,a') = (1 — aa')/2 
so that there is precisely one hydrogen on each lattice edge. The model 
considered by Attard is described by 

E{a, a') = w+, a = a = 1 
= W-, (T = a' = — 1 

= 1, a = -a'. (4) 

The percolation model of ^ is equivalent to a special case of (§) with = 
w_ = e"^^ and 

E{a,a') = 6^^'^'^'+^) = e^(coshir)(l + W), (5) 

where z = tanhK. For our purposes, we shall restrict our considerations to 
the percolation model (|]). 
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Figure 2: The 8-vertex model configurations. 

Attard and Batchelor P, |^ adopted an arrow representation for tlie liydro- 
gen configurations which, due to the occurrence of defects, led to a 14- vertex 
model with Bjerrum defects. A weak-graph transformation is then carried 
out for the 14- vertex model. Here, we expand the partition function directly. 
Substituting (|) into (|^) and expanding the second product over the edges 
of £, we obtain an expansion of 2^^ terms. To each term in the expansion 
we associate a bond graph by drawing bonds on those edges corresponding 
to the z factors contained in the term. This leads to a 16- vertex model on L. 
Besides an overall Boltzman factor (e^ cosh i^)^^, the 16- vertex model has 
vertex weights 

where the product is taken over those incident edges with bonds. The sym- 
metry relation (|]) now implies that W = ^ whenever there are an odd number 
of incident bonds, and the 16-vertex model becomes an 8-vertex model. 

Using the bond configurations of the 8-vertex model shown in Fig. 2, it 
is straightforward to deduce using (||) the following vertex weights: 

Wx = u {ai,a2, (73, a4,) = 2{a + b + c) 

W2 = 2:^ ^ai(T2(T3Cr4Co'(cri, (72, 0-3, (T4) = 22:^(0 + 6 + c) 

W3 = Wi = Z^a2a4,uj{ai,a2,a3,a4) =2z{-a-b + c) 
W5 = We = z'^aia2Uj{ai,cr2,cr3,a4) =2z{-a + b-c) 
W7 = Ws = z'^a2(T3Uj{ai,(T2,(T3,cr^) = 2z{a - b - c), (7) 

where Wi is the vertex weight of the ith vertex and we have used the fact that 
in the non-vanishing tu's we have (7i(T20"3(T4 = 1. Now, in an 8-vertex model 
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configurations, vertices 3 and 4, 5 and 6, and 7 and 8, always occur in pairs 
and/or in even numbers. Therefore we can conveniently replace the relevant 
weights by their absolute values and arrive at, after dividing all weights by 
a common factor 



Wi = 




(a + b + c) 


W2 = 


z{a 


+ b + c) 


Ws = 


W4 


= \— a — b + c\ 


w, = 


We 


= \— a + b — c\ 


Wj = 


Ws 


= \a — b — c\. 



The vertex weights describes an 8-vertex model in an external electric 
field h = {\nz)/2 in both the vertical and horizontal directions fl]]. 

More generally, if we allow different values of w+ = W- = Wi,i = 1,2 
in (^) and for the horizontal and vertical edges respectively, and write 
h = (lnzi)/2 and v = (ln2;2)/2, where Zi = {wi — l)/{wi + = 1,2, then 
one arrives at an 8-vertex model with weights 



Wi 


= e^-^^ia + b + c) 




= e-''~''{a + b + c) 




= e^-^l-a-b + c] 




= e'^-^l - a - b + c\ 


w^ 


= Wq = \ - a + b - c\ 


W-j 


= Ws = \a - b - c\. 



In ensuing discussions we shall consider the general model (9). 



3 The free-fermion solution 

The free-fermion model is defined as a particular case of the eight-vertex 
model in which the vertex weights satisfy the relation 

W1W2 + W^Wi = W^Wq + VTrWg, (10) 

a condition equivalent to the consideration of a noninteracting many-fermion 
system in an 5- matrix formulation of the 8-vertex model [|10|. In the present 
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case the free-fermion condition (10) is satisfied when either a = or 6 = . 
Without the loss of generality, we consider a = and b > c. 

The closed expression for the free energy of the free-fermion model is 
well-known Q and after some algebraic manipulation using results of [^], we 
obtain 

-/3/= hm llnZ = ln(6 + c) + -^ r^ci01n(A + gi/2), (11) 
N^oo i\ 47r Jo 

where 

Q = [sinh(2t; + 2/;,) + fc^sinh (2t; - 2/i) + 2A;cosh2t;cos</)]^ + 4fc^sin(/)^, 
A = cosh {2v + 2h) +P cosh {2v ~2h)+2k sinh 2t; cos (12) 

with k = {b — c){b + c). 

The critical condition of the free-fermion model is given by 

Wi + W2 + W^ + Wi = 2 max {Wi, W2, W^s, W^}, (13) 

where 





= e'^+''(6 + c) 


W2 


= e-^-^b + c 




= e^-\b~c) 




= e^'^ib-c). 



(14) 



Thus, the h-v plane is divided into four regions depending on which vertex 
1, 2, 3 or 4 has the largest weight. Denoting the four regions by /, //, /// 
and IV respectively, as shown in Fig. 3, the critical condition (13) can be 
rewritten as 



b 


tanh V + e^^ 


c 

b 


1 — e^^ tanh v 
1 — e^^ tanh v 


c 

b 


tanh V + e^^ 
tanh V — e^^ 


c 

b 


1 + e^^ tanhv 
1 + e^^'tanhu 


c 


tanh V — e^^ 



region /, 

region //, 

region ///, 

region IV. (15) 
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2d=ln[(b+c)/(b-c)] 
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Figure 3: The four regions in the plane h-v denoted by /, //, /// and IV 
are, respectively, the regions where vertex 1, 2, 3 or 4 has lowest energy. 

The critical condition (15) is plotted in Fig. 4a — Ad for four specific values 
of b/c. Generally, the free energy exhibits a logarithmic singularity at the 
phase boundaries with exponents a = a' = . When b = c for which Q is 
a complete square, however, we have 

-(3f = ma.x{\nWu\nW2} 

= \n{2b) + \h + v\ (16) 

and the phase boundary h + v = separates the two frozen states Wi and 

4 The Bethe-like lattice 

Now, we consider the 6- vertex model with bond defects on a Bethe-like lattice 
with plaquettes as shown in Fig. 5. The study of systems in the Bethe-like 
lattices is an alternative approach to the usual mean-field theory. The main 
features of the model under investigation will be obtained by studying the 
properties of the free energy. 

A free energy in a region deep inside a Bethe-like lattice must be carefully 
defined. It cannot be obtained by directly evaluating the logarithm of the 
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Figure 4: Phase boundaries of the free-fermion modeL 



partition function in which the contribution from the outside of this region 
is not neghgible and as result the system exhibits an unusual type of phase 
transitions without long-range order |]n], |TB|. Recently, a method for the 
surface independent free energy calculation is presented |15|. The free 
energy /□ per plaquette of our model is expressed as 



- (3fn = lim ^ (In Z„ - 3 In 



(17) 



where Zn and Zn-i is the partition functions of the 6- vertex model with bond 
defect on the Bethe-like lattice consist of n and n — 1 generations respectively. 

The calculation on the Bethe-like lattice is based on a recursion method 



1^. When the tree is cut at the central plaquette, it is separate into 4 
branches, each of which contains 3 branches. Then the partition function of 
interest (|3f) can be written as follows: 



N, 



(n) 



^(^01, O-02, (T03, O-04)5'n(o"0l)fl'n(o-02)fl'n(o-03)5'n(o-04), 

Woi} 



(18) 
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Figure 5: The Bethe-like lattice. 



where iV^^"^ = 2(3" — 1) is the number of bonds, n is the number of generations 
and Qni^'oi) is in fact the partition function of a branch nearest to the Oi site. 

Each branch, in turn, can be cut along any site of the first generation, 
then the expressions for Qnio'oi) can therefore be written in the form: 

9n+i{croi) = ^{1 + 2;croiO"i3)^^(o-ii, 0-12, cri3, ai4,)gn{crn)gnicri2)9n{cru) (19) 
Wii} 

After dividing gn{—) by gn{+), we obtain a recursion relation for Xn = 
gn{—)/gn{+)- Let us consider the case when series of solution of recursion 
relation converge to a stable point at n — *• oo, namely, 

lim Xn = X. 

n— >oo 

We obtain the following equation: 

(1 - z)x'^ + (1 +z)x 



X 



[1 + z)x'^ + (1 - z)x' 



(20) 



We are now in a position to compute the free energy per plaquette of our 
model. Using Eqs. (0), ([T^), (|T^ and (^), the expression for the free 
energy functional can be written as 



lim ifiV^") 



2 

1 r 



(n-l) 



In 



UJ+l 



(21) 



+ lim 

n^oo 2 



In $(x(„_i)) + In ^(x(„)) - 3 In ^(x(„_i)^ 
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where 

= 2(a + 6 + c)x^ and 



It easy to see that A*", 



in) 



3N, 



(n-l) 



--{a + b + c) + z)x + l- zf (22) 
4 for all n. Thus the free energy per 



plaquette can be finally written as 



-5/, = ,„(i;i±i)!(a + t + c) + l„(i±i^^±i^. (23) 

Together with the expression (^) for x it gives the free energy per plaquette 
of the 6-vertex model with bond defects. 

The equation of state (|20|) always has x = 1 as fixed point solution. In 
this case the free energy per plaquette is 



Pfa = \n ^ Z ' {a + h + c) 



(24) 



In the case ofa = 6 = c=l, we recover the result obtained by Attard 
and Batchelor P] for the six-vertex model with bond defect in the mean-field 
(independent vertex) framework, 

-/3/n=ln^(w; + l)2, 



which reduces to Pauling's estimate when w = p 
For the specific heat, we obtain 



c = A 

dT 



dT 



2w 



or, explicitly, 

1 2^2 



{\nwf + 



ahi\n.a/hY + ac(lna/c)^ + hciXnh/c) 



C 



T2 (coshis:/T)2 
where Eab = ea-eb, e 



+ 



(a + 6 + c)2 
T2(l + e'^"/^ + e^-/^)2 



(25) 



(26) 



ac — ~ and 



-ea/T 



-Sc/T 



= 6""'", 6 = 6"'"", c = e 
The specific heat versus T /K for a = 6 is plotted in Fig. 6. 
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Figure 6: The specific heat for the Bethe-hke lattice of Fig. 5 as function 

of T/K for a = 6 and for p={ea- e^jK = 0, 1. 

5 Summary and discussion 

In this paper we have considered a 6- vertex model of hydrogen-bonded crys- 
tals with bond defects. We have established the exact equivalence of the 
6-vertex model with bond defects with an 8-vertex model in an external 
electric field. Using this equivalence we solve exactly our model in the free- 
fermion subspace. We also obtain the exact solution of the 6-vertex model 
with bond defects on a Bethe-like lattice. 

In P, one of us has used the percolation representation of the hydrogen- 
bonded model to argue that the specific heat of the system will increase as the 
temperature decreases from very high temperature. Figure 6 indeed shows 
such behavior. However, the specific heat in Fig. 6 does not diverge. It is of 
interest to calculate the specific heat for the square or diamond lattices by 
Monte Carlo method. We expect to find singular behavior of specific heat in 
such systems. 
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